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Sununary 

This  article  examines  the  dynamic  plastic  response  of  an 
impulsively  loaded  ideal  fibre-reinforced  (strongly  anisotropic) 
rigid-plastic  beam.  It  is  demonstrated  that  the  response  duration 
and  permanent  transverse  displacements  are  significantly  less  than 
the  corresponding  quantities  in  an  "equivalent"  rigid  perfectly 
plastic  isotropic  beam. 
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Notation 


x,y,z, 


location  of  plastic  hinge  (Figure  2(b)) 

(Q1/m) 1/2 

mass  per  unit  length 

strain  hardening  index  in  equation  (1) 

V°o 

time 

response  duration 

response  duration  for  a rigid  perfectly  plastic 
isotropic  beam 
displacement  in  z sense 
permanent  displacement  in  -z  sense 

maximum  permanent  displacement  in  a rigid  perfectly 
plastic  isotropic  beam 

cartesian  coordinates  defined  in  Figure  1 
beam  thickness 


beam  span 


V°1 


0oH  / 4 

transverse  shear  force  defined  in  Figure  1 
parameters  defined  by  equation  (1) 
transverse  shear  force  defined  by  equation  (1) 
impulsive  velocity 


shear  strain 


VD/C 

uniaxial  yield  stress 

JL  ( ) 

at  ' ’ * 
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Introduction 

t 

The  dynamic  plastic  response  of  some  ideal  fibre-reinforced 
(strongly  anisotropic)  rigid-plastic  beams  were  examined  recently 
in  References  [1]  to  [6]  using  a continuum  theory  which  made  no 
distinction  between  the  behavior  of  the  fibres  and  the  response  of 
the  matrix.  These  theoretical  predictions  were  particularly 
attractive  because  the  main  features  of  response  were  obtained 
with  little  effort.  Moreover,  it  was  demonstrated  in  Reference 
[2]  that  the  response  durations  and  permanent  transverse  displace- 
ments of  several  ideal  fibre- reinforced  beams  loaded  dynamically 
were  less  than  the  corresponding  values  for  similar  rigid-plastic 
isotropic  beams. 

The  beams  examined  in  all  the  previous  References  ([l]-[6]) 
were  impacted  transversely  by  a point  mass  M*  travelling  with  an 
initial  velocity  V . It  is  the  purpose  of  this  present  note  to 
seek  the  dynamic  plastic  response  of  a beam  which  is  subjected  to 
an  impulsive  transverse  velocity  VQ  distributed  uniformly  over  the 
entire  span.  The  particular  case  of  a clamped  beam  is  studied 
herein  because  the  experimental  arrangement  required  for  tests  is 
fairly  straightforward  and  comparisons  may  be  made  with  the 
corresponding  theoretical  and  experimental  results  which  have  been 
reported  in  the  literature  for  the  dynamic  plastic  response  of 
similar  isotropic  beams  [7]. 

+ The  point  mass  M is  infinite  in  References  [1]  and  [3]  as 
discussed  in  Reference  [2]. 


Basic  Equations 


The  theoretical  foundations  on  which  the  analysis  in  this 
note  rests  are  described  in  References  [1]  and  [2]  to  which  a 
reader  is  referred  for  further  details.  The  material  is  assumed 
to  be  incompressible  and  inextensible  in  the  direction  of  the 
fibres  along  a beam  axis.  In  addition,  the  deformations  are 
assumed  to  remain  infinitesimal,  while  material  elasticity  and 
transverse  wave  propagation  effects  are  neglected:  simplifications 
which  are  customarily  made  in  theoretical  work  on  the  dynamic 
plastic  response  of  isotropic  beams  [8-10].  These  assumptions 
imply  that  only  shear  deformations  can  develop  in  ideal  fibre- 
reinforced  beams. 

The  material  is  assumed  to  be  rigid  strain  hardening  [3]. 

Qp  = ±{Qq  + Qx ! Y I n>  (1)+ 

where  y = 3w/3x  is  the  shear  strain  and  Qq,  and  n are  positive 

constants.  Plastic  flow  of  the  material  is  controlled  by  the 
magnitude  of  the  transverse  shear  force  Qp  as  discussed  in 
Reference  [2],  and  the  material  remains  rigid  during  unloading. 

The  equation  of  transverse  equilibrium  for  a beam  is 

3Q/3x  ♦ P(x,t)  - m32w/3t2  (2) 

where  the  various  symbols  are  defined  in  Figure  1 and  in  the 
accompanying  notation.  If  an  element  of  a beam  at  x = a is 
crossed  by  a velocity  discontinuity  which  travels  with  a velocity  a then 
the  conservation  of  momentum  (dynamic  jump  condition)  demands 

[Q(a#  t)  ] - - ma[3w/3t]  (3) 

+ The  theoretical  analyses  in  References  [1]  and  [2]  were 

developed  for  a rigid  linear  strain  hardening  material  with 
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at  x = a,  while  continuity  of  displacement  (kinematic  jump 
condition)  requires 

a[3w/3x]  = - [3w/3t]  (4) 

at  x = a,  where  [x]  means  X(a+)  - X(a  ) and  denotes  the  difference 
in  X on  either  side  of  the  discontinuity  at  x = a(t).  In  passing/ 
it  should  be  remarked  that  equations  (3)  and  (4)  are  identical  to 
equations  (19a)  and  (21a)  of  Symonds  [11]. 

Clamped  Beam  Loaded  Impulsively 

Consider  a stationary  ideal  fibre-reinforced  rigid-plastic 

beam  which  is  subjected  to  time  t = 0 to  a velocity  VQ  distributed 

uniformly  over  the  entire  span  2L  as  shown  in  Figure  2(a).  Large 

transverse  shear  forces  (Q  ) will  develop  immediately  at  the  fully 

P 

clamped  supports  and  give  rise  to  a slope  discontinuity  (hereafter 
called  a "shear  hinge")  according  to  equation  (1) . These  shear 
hinges  propagate  with  a velocity  a inwards  from  both  supports 
towards  the  beam  centre  as  indicated  in  Figure  2 (b) . Plastic 
deformation  is  confined  to  the  two  travelling  shear  hinges,  while 
the  central  portion  of  the  beam  -a  £ x £ a moves  as  a rigid  body 
with  a velocity  V(t)  and  the  two  outer  regions  a <_  x £ L and 
-L  < x £ -a  remain  stationary.  This  problem  is  symmetric  about 
x = 0 so  that  only  the  region  0 < x < L is  considered  henceforth. 
Now,  equation  (3)  becomes 

ma[ 3w/3t]  = -[Qp]  (5) 

across  a travelling  shear  hinge , which  using  equation  (1)  can  be 


written 


ma[ 3w/3t]  * + Qx  [ | Y | n] . 


(6) 
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Equations  (4)  and  (6)  give 

a2  = ± c2 [ | y | n]/[ Y 3 • 


(7) 


1/2 


where  a is  the  velocity  of  a shear  hinge  and  c = (Q^/m) 
is  the  propagation  speed  of  a shear  hinge  in  a linear  work 
hardening  material  (n  = 1) . It  is  evident  that  the  velocity  of 
propagation  a is  not  constant  when  n / 1 since  y is,  generally 
speaking,  a function  of  x. 

If  the  equation  of  motion  (equation  (2)  with  P = 0)  in  the 
central  region  0 < x < a is  integrated,  then 


maV  = -Q 


(8) 


since  Q = 0 at  x = 0 and  Q = Q at  x = a,  where  a refers  to  the 

o 

left  hand  side  of  the  shear  hinge  at  x = a in  Figure  2(b). 
Furthermore,  equations  (4)  and  (6)  give 


V = -ay (a) 


(9) 


and 


maV  = -Qi | y (a+) |n  (10) 

respectively,  where  a+  refers  to  the  right  hand  side  of  the  shear 
hinge  at  x = a in  Figure  2(b).  Equations  (9)  and  (10)  may  be 
recast  in  the  form  ^ 2 

Y (a+)  - (m/Qi)n*1vn'fl  (11) 

and  _L_  n^i 

• n 4- 1 n4l 

(12) 


a = - (Qi/m)n+1 V n+1 


which  also  may  be  obtained  from  equations  (7)  and  (11).  Now, 
equation  (8)  can  be  written 


m 


{L  + 


r* 


adt  }_1 


V 


(13) 
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or 


n-1 


-Q„ 


V = 


m 


■ Q •>  n+1  r t n+1 

{L-hn  J v(t)  dU 


-1 


(14) 


when  using  equation  (12) 

If  equation  (12)  is  divided  by  equation  (8) , then 

2 n-1 


da  _ ma  c 
dv  - ^ 


n+1  n+1 

V 


(15) 


which  has  the  solution 


_ _ T r (n+1) m c 

a - L exP*  {"~2nQT_ 


2 

n+] 


2n 

n+1 


2n 

n+1 


-V0 


(16) 


since  V = VQ  and  a = L when  t = 0.  The  final  position  of  the 
travelling  shear  hinge  (a^)  is  given  by  equation  (16)  with 


V = 0 , i.e.  , 


a,  = L exp. { 


- (n+l)m  c 


2 2n_ 

n+1  n+1 


‘f  “ 1 2nQ  o 

o 


} 


(17) 


Now,  substituting  V from  equation  (16)  into  equation  (11)  gives 

2nQ . 


2n  h 


Y (a+)  = {- 


(n+1)  Q 


r 1 V n 

, a o n+1 

loge  L + l— J } 


(18) 


which  may  also  be  interpreted  as  the  permanent  shear  strain 

since  all  plastic  deformation  is  confined  to  the  travelling  shear 

hinges.  Thus,  the  permanent  transverse  displacements  in  the 

L 

region  af  < x <L  are  wf(x)  = j y(x)  dx++,  or 


It  is  evident  from  equation  (16) that  the  travelling  shear 
hinge  never  reaches  the  beam  centre. 

++wf(x)  is  defined  as  positive  in  -z  sense. 


d 
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since  V decreases  as  the  shear  hinge  travels  inwards  towards  the 
beam  centre  according  to  equation  (16) . It  is  apparent  that  the 
foregoing  theoretical  predictions  are  both  kinematically  and 
statically  admissible  and  therefore  exact  within  the  spirit  of 
the  theory. 

The  transverse  shear  force  in  the  outer  zone  a < x < L is 

Q = Qp(y(a+))  and  therefore  the  associated  bending  moment 

f * + 

distribution  according  to  M = Q dx  is  M = M - Q (y  (a  ) ) (L-x) , 

* 

where  M is  the  magnitude  of  the  bending  moment  at  the  supports. 
Similarly,  Q = QQx/a  and  M = M*1  + QQ(x2-a2)/2a  in  the  inner 
region  o£  x < a,  where  M*1  is  the  value  of  the  bending  moment 
at  the  travelling  plastic  shear  hinge.  However,  the  bending 
moment  must  be  continuous  across  a travelling  shear  hinge  when 
rotatory  inertia  is  neglected  [11].  Thus,  M*1  = M - (y  (a+)  ) (L-a) 
Now,  the  bending  moment  (M)  does  not  enter  into  the  yield 
condition  because  an  ideal  fibre-reinforced  beam  is  inextensible 
along  its  length  and  therefore  any  value  of  M required  for 
equilibrium  may  be  supported  without  any  associated  strains  and 
without  penetrating  the  yield  surface.  It  is  evident,  therefore, 

that  the  foregoing  theoretical  solution  is  valid  for  any  value  of 

* 

M at  the  supports.  In  particular,  the  theoretical  predictions 

* 

are  valid  for  "simple"  supports  with  M = 0 as  well  as  the  fully 
clamped  case. 


Impulsive  Loading  of  a Clamped  Beam  with  Linear  Strain  Hardening  (n=l) 
The  theoretical  analysis  in  the  previous  section  is  particularly 
simple  when  n = 1 since  the  shear  hinge  then  propagates  with  a 
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constant  velocity  a = -c  according  to  equation  (12) . 

Integration  of  equation  (13)  with  n = 1 and  V = V when 

o 

t = 0 gives 

V - Vo  + si  l°9e<l  - f > (22) 

which  predicts  a duration  of  response 

T T -mcV  1 

fcf  = c L1  ‘ eXp*  { ~Q^°  }J  (23) 

when  V = 0.  Thus,  the  final  position  of  the  shear  hinge  is 
af  = L - ctf , or 

- me  V 

af  = L exp. { — } . (24) 

° L 

The  permanently  deformed  profile  of  the  beam  is  wf  = I yf(x)dx+ 

•*x 

which  when  using  yf(x)  = y(a+),  equation  (11)  with  n = 1 and 
equation  (22)  gives 

"f  - ¥ (1  - E»  - ^ + - E > (25> 

1 

for  af  £ x < L and 
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Discussion 

It  is  evident  from  equation  (29)  that  the  final  position  of 

a shear  hinge  (a^/L)  is  near  the  supports  when  B = Xq  is  small 

and  n = 1.  In  this  circumstance,  a large  central  portion  of  a 

beam  remains  undeformed  and  the  plastic  deformation  is  confined 

to  the  vicinity  of  the  supports  as  shown  in  Figure  3.  However, 

the  size  of  the  deformed  zone  near  the  supports  of  a beam  increases 

as  the  strain-hardening  index  n is  decreased,  while  the  associated 

permanent  transverse  deflections  decrease  according  to  equation 
% 

(19)  and  as  indicated  in  Figure  3.  The  influence  of  the  material 
strain  hardening  parameter  (q)  and  strain  hardening  index  (n)  on 

the  maximum  permanent  transverse  displacements  according  to 
equation  (19)  is  shown  in  Figure  4,  where  the  abscissa  is 
X (q) 1/2  = VQ(m/Qo)1/2. 

The  permanent  transverse  deflection  at  the  mid-span  of  a fully 
clamped  rigid  perfectly  plastic  isotropic  beam  loaded  impulsively 
with  a uniform  velocity  VQ  is  [12] 

2 2 

wfi  - mVQ  L /(6Mq)  (32) 

2/ 

where  MQ  ■ oqH  4 and  H is  the  thickness  of  a beam  with  a rectangular 
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cross-section.  If  plastic  yielding  of  a cross-section  is 

controlled  by  the  Tresca  yield  condition,  then  the  fully  plastic 

transverse  shear  force  is  Q = a H/2  and  therefore  M = Q H/2. 

o o o o 

Thus,  equation  (32)  may  be  written 

2 

wfi/L  = A q L / ( 3H)  (33) 

which  when  divided  by  equation  (30)  for  n = 1 gives 


§ (L} 

3(6  - 1 + e“B)  H 


where  6 = Aq.  Equation  (34)  predicts 

wfi/wf  = 2L/(3H) 


(35) 


when  8 <<  1 (e.g.,  very  small  strain  hardening  parameter  Q^) . 
the  variation  of  wf^/wf  L/H  and  6 according  to  equations 

(34)  and  (35)  is  shown  in  Figure  5.  These  results  indicate  that 
the  permanent  transverse  displacements  of  this  ideal  fibre- 
reinforced  (strongly  anisotropic)  rigid  plastic  beam  with  n = 1 


are  much  smaller  than  those  of  an  "equivalent"  rigid  perfectly 
plastic  isotropic  beam,  even  when  ■*  o (i.e.,  6 -*■  o)  . 

The  response  duration  of  a fully  clamped  rif'id  perfectly 
plastic  isotropic  beam  is 

fcfi  " mVoL  /(4Mo)  ' (36) 

while  equation  (23)  or  equation  (28)  predict 

tf  = m V L/Q  (37) 

t o o 

when  8 <<  1 (e.g.,  very  small  strain  hardening  parameter  . 

Thus , 

tfi/tf  = V(2H)  (38) 


i 
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which  indicates  that  the  response  duration  of  an  ideal  fibre- 
reinforced  (strongly  anisotropic)  rigid-plastic  beam  is 
considerably  less  than  an  "equivalent"  rigid  perfectly  plastic 
isotropic  beam. 

Conclusions 

This  article  examines  the  dynamic  plastic  response  of  an 
impulsively  loaded  ideal  fibre-reinforced  (strongly  anisotropic) 
rigid-plastic  beam  with  material  strain  hardening.  It  is  demonstrated 

that  the  response  duration  and  permanent  transverse  displacements 
are  significantly  less  than  the  corresponding  quantities  in  an 
"equivalent"  rigid  perfectly  plastic  isotropic  beam.  Thus,  it 
appears  that  a potential  for  considerable  weight  savings  may  exist 
for  energy  absorbing  systems  made  from  materials  characterised  as 
ideal  fibre-reinforced  rigid-plastic.  However,  experimental 
investigations  are  required  to  establish  if  the  ideal  fibre- 
reinforced  material  model  is  useful  for  strongly  anisotropic 
beams  loaded  dynamically.  The  particular  case  examined  herein 
offers  an  attractive  experimental  arrangement  for  experimental 
work  since  similar  rigid-plastic  beams  have  been  investigated  by 
many  authors. 
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Figure  3. 


Sign  Convention. 

(a)  Impulsive  Loading  of  a Fully  Clamped  Beam 

(b)  Transverse  Displacement  Profile  at  time  t. 
Permanent  Transverse  Displacement  Profiles  according 
to  Equation  (19)  with  X = 0.1  and  q = 4. 
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® : n = 1,  af/L  = 0.6703,  ct./L  = 0.3297,  w /L  = 0.01758 
^ r max 

n = 1/3,  af/L  = 0.07967,  ctf/L  = 0 . 00635 , wmax/L  = 0 . 01259 

(D  : n = 1/5,  a./L  = 0.00381,  ct./L  = 0.00152,  w /L  = 0.01055 
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Figure  4.  Influence  of  Material  Strain  Hardening  Parameter 
(q  = Qj/Q  ) and  strain  hardening  index  n on  the 
Maximum  Permanent  Transverse  Displacements  of  Impulsively 
Loaded  Beams  According  to  Equation  (19) . 

Figure  5.  Comparison  of  Maximum  Permanent  Transverse  Displacements 
of  Rigid  Perfectly  Plastic  Isotropic  (wfi)  and  Ideal 
Fibre-Reinforced  (Strongly  Anisotropic)  Rigid  Plastic 
(wf)  Beams  Loaded  Impulsively  According  to  Equations 
(34)  and  (35). 
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